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ABSTRACT 

Aims. In this paper, we study the properties of solar granulation in a facular region from the photosphere up to the 
lower chromosphere. Our aim is to investigate the dependence of granular structure on magnetic field strength. 
Methods. We use observations obtained at the German Vacuum Tower Telescope (Observatorio del Teide, Tenerife) using 
two different instruments: Triple Etalon SOlar Spectrometer (TESOS), in the Ban 4554 A line to measure velocity and 
intensity variations along the photosphere; and, simultaneously, Tenerife Infrared Polarimeter (TIP-II), in the Fel 1.56 
/xm lines to the measure Stokes parameters and the magnetic field strength at the lower photosphere. 
Results. We obtain that the convective velocities of granules in the facular area decrease with magnetic field while the 
convective velocities of intergranular lanes increase with the field strength. Similar to the quiet areas, there is a contrast 
and velocity sign reversal taking place in the middle photosphere. The reversal heights depend on the magnetic field 
strength and are, on average, about 100 km higher than in the quiet regions. The correlation between convective velocity 
and intensity decreases with magnetic field at the bottom photosphere, but increases in the upper photosphere. The 
contrast of intergranular lanes observed close to the disc center is almost independent of the magnetic field strength. 
Conclusions. The strong magnetic field of facular area seems to stabilize the convection and to promote more effective 
energy transfer in the upper layers of the solar atmosphere, since the convective elements reach larger heights. 
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1. Introduction 

Solar faculae are areas on the solar surface surround- 
ing active regions and appearing bright toward the limb. 
Understanding their properties is important from several 
points of view. On the one hand, facular contrast influ- 
ences total solar irradiance variations. On the other hand, 
faculae are believed to be composed of conglomerates of 
magnetic elements. Therefore, in faculae we can observe 
the interaction of convection with relatively strong mag- 
netic fields, which is interesting in itself from a physical 
point of view and also to constrain state of art numerical 
models of magneto-convection. 

It is well established that the origin of facular bright- 
ness excess is related to the presence of magnetic field: the 
brightness excess is caused by the modified radiative trans- 
fer throu gh magnetic atmospheres viewed from different an- 
gles f e.g. iKeller et al.L [20 04: Carl sson et all . 120041: IVogleil 
l2005H Ok unev fc Kneeii 12005; .Steijr2005[ ). The variation 
of facular brightness near the solar limb can provide in- 
formation about the structure and properties of the mag- 
netic elements composing facular regions. Measurements 
of centre-to limb variations of the contrast done by dif- 
ferent authors mos tly agr e e that the contrast increases 
till u = 0.2 - 0.4 (iMulled. [19771 lAuffret fc Mulle'R [1991 



ToDka ct al.', '1997*: Orti z et all. 120021: lOkunev fc Kneer 



2004: .Hirzbcr gcr fc Wichrl. 120051 ). A controversy exists for 



the extreme li mb, where it is not clear if the contrast fur - 
ther increases (|de Boer et all . 119971: ISiitterlin et"aI1 . 119990 



or falls off toward the limb. The variety of results can 
be attributed to the difference in spatial resolution of 
observations, facular size, its magnetic field strength, as 
well as selection criteria. At high spatial resolution (0."1- 
0."2) and at disk centre, facular regions appear as con- 
glomerates of brigh t points, small pores and smal l -scale 
granular structures (iLites et all . 120041: iBerger et all . 120071: 
llMaravan fc Scharmerl . l2010l ). The contrast of faculae at disk 
centre is slightly positive or close to z ero, though some 
negative values ar e also detected (see iTitle et all . 119921 : 
iTopka et al.L Il997l and references therein). 

Several classes of models claim to reproduce facular 
brightnes s prope r ties. T he "hot wall" model, initially pro- 
posed bv ISpruifl (|1976[ ). consists of a magnetic flux tube, 
evacuated due to the presence of magnetic field, similar 
to Wilson depression in sunspots. Viewed from an angle, 
the linc-of-sight penetrates deeper into the tube because 
of its lower density and, thus, hotter layers become visi- 
ble. T h is idea was fo ll owed in works by e.g., Topka et al.l 
(119970 : ISteine^ (l2OO50 :[()k unev fc Kneed doO.C " who pro- 
vide evidences that the "hot wall" model can closely re- 
produce observed facular properties, such as continuum 
contrast, its centre-to- limb variation, as well as the depen- 
dence of contrast on the magnetic field. Much more complex 
models o f facular reg i ons based on "realist i c" MHD sim- 
ulatio ns (|Keller et al.l . 120041: ICarlsson et all . 120041: IVogled. 
|2005| ) suggest that faculae are seen bright on the limb be- 
cause hot granular walls become visible through transpar- 
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ent magnetic field concentrations. The three-dimensional 
structure of granules becomes apparent at faculae near the 
solar limb in high resolution observations and is well repro- 
duced in the simulations. Two factors are crucial for the 
facula formation: the shape of the granule limb- ward of the 
flux concentration and the size of the magnetic field concen- 
tration. If the flux concentration is too small, the opacity 
reduction is not sufficient to have the continuum intensity 
formed exclusively in the hot granule. A quantitative dis- 
agreement still remains in the values of the peak brightness 
in faculae, being subs tantially larger in simulations com- 
pa red to observations (jKeller et al.l . |2004|) . 

iBerger et al.l (|2007l ) question the dependence of granu- 
lar bri ghtness on magnetic field strength obta i ned in earlier 
works (jTopka et al.l . ll997l: lOrtiz et all l2002[ ). iBerger et all 
(|2007f) analyze extremely high spatial resolution data (close 
to 0."1) obtained at the Swedish Solar Telescope on 
La Palma. They claim that if, instead of analyzing the 
brightness for binn ed magnetogram signals (as is done in 
iTopka et al.l Il997[ ) , one analyzes magnetic fiux density for 
facular points segmented from the data set, there is no de- 
pendence of t he br ightness on the magnetic fiux density. 
IBerger et al] (|2007t ) propose that what we see as faculae 
are granular walls, and not interiors of the magnetic fiux 
tubes. Since granules all have similar properties, there is no 
dependence of their brightness on the magnetic field, and 
the magnetic field only plays an indirect role, making the 
atmosphere transparent in front of the granules. 

The strong magnetic field of facular regions modifies 
the properties of convection. Simulations of magnetic fiux 
emergence through the convection zone show that dur- 
ing the emergence granules become la rger, more elon- 
gated and smoother (ICheung et al.l . [2008f ). Observationally, 
in already formed facular region granulation presents a 
lot of fine structuring at high spatial resolution: isolated 
bright points, strings of bright points and dark micro-pores , 
ribbo n s, or more c ir cular fiower structures (iTitle et aLl . 
Il992t IBerger et aP . 12004 iNaravan fc Scharmeij |2010D . 
Gran u les become smaller than in quie t areas (jMulleil . 
[1971 iSchmidt et all . [1981 iTitle et all . [1991 and in- 
tergranular lanes are characterized by the presence of 
micro-pores. Granular velocities in facular and plage 
areas are found to be generally lower than in the 
nearby quiet areas (iNesis fc Mattid . 119891: iTitle et al.l . 119921: 
INaravan fc Scharmeil I2010D . Velocities measured in plage 
and facular areas are found to depend on the magnetic 
field strength. From medium-to-high (430 km) resolu- 
tion ob servations o f the disk centre plage in Nil 6768 
A line, ITitle et al.l (jl992[) obtain an increase of down- 
flow velocities with magnetic field strength up to 600 
G, and a decrease for stronger fiel d s. A similar result i s 
reported by JX lont agiTe et al.l (119961): IBerger et al.l i2004 ; 
iMorinaga et al.l (|2008[ ): INaravan fc Scharmeil (|2010l ). The 
correlation between velocity and intensity fields, typical for 
granulation, is found t o be partially destroyed by the mag - 
netic field jRimmeH, 12004 INaravan fc Scharmeil |2010D . 
O ne of the po ssible reasons of t he de-correlation, proposed 
by INaravan fc Scharmeij (|2010[ ) is that small scale convec- 
tion in plage areas does not overshoot to the same height 
as field-free convection and leaves only weak traces at the 
line-forming region. 

It is clear that further observational studies of magneto- 
convection in strongly magnetized plage and facular ar- 
eas are needed to constrain the relationship between 



magnetic field and granulation, as well as the proper- 
ties of magnetic elements. Here we report on such an 
observational study based on state-of-the-art simultane- 
ous observations done at the German Vacuum Tower 
Telescope at the Obser vatorio del Teide (T enerife) with 
two instruments: TIP-II (jCollados et al.l . l2007f ) and TESOS 
([Tritschler et al.l . l2002l) . We analyze correlation between 
granular /intergranular intensities, velocities, and magnetic 
field; as well as the results of the heights of sign reversals 
of contrast and velocity as a function of the magnetic field. 



2. Observations 

The observations were obtained on the 13th of November, 
2007, at the German Vacuum Tower Telescope at the 
Observatorio del Teide (Tenerife). We simultaneously ob- 
served three spectral regions: Fei 15643-15658 A, Ban 4554 
A, and Can K 3968 A. 

Using the filtergrams in Can K, we have selected a fac- 
ular area located close to disk centre. There was no sunspot 
activity at the time of observations and the observed facular 
area was not attached to any active region. 

The acquisition of the data was co ntrolled by the 
Tenerife Infrared Polarimeter (TIP-II, see ICoUados et al.l 
l2007l ). to scan repeatedly a small area of the solar surface. 
The slit length was 84" and the slit width 0."5. The pixel 
size of the CCD of TIP was 0."185. One scan consisted 
of 15 steps, separated 0."35, with a total scanned area of 
84"x5."5. The time interval between successive scan steps 
was 27.3 sec. Thus, it took 6 min 50 sec to record each 
complete series of 15 steps. All four Stokes parameters were 
recorded in the 1.56 ^m region. In total, 22 scans were done. 
The spectral sampling of the infrared data was of 14.7 mA. 

Simultaneously with TIP-II, the TESOS instrument 
was making scans through the spectral profile of the Ban 
4554 A line. TESOS (Triple Etalon SOlar Spectrometer; 
iTritschler et"all . |2002[) is a Fabry-Perot interferometer 
equipped with three etalons which allow to scan spectral 
lines in a two-dimensional field of view. The CCD of TESOS 
had 512x512 pixels of 0."089 size. Thus, the total size of 
the observable area on the Sun was 45."6x45."6. The Ban 
4554 A line was covered by 37 spectral points with a spec- 
tral sampling of 16 mA. During the scan, TIP was moving 
the whole field of view of TESOS at steps of 0."35. However, 
at every moment of the observations, the area scanned by 
TIP and the field of view of TESOS overlapped. Each com- 
plete scan of the 37 spectral points of the Ba n line profile 
took 25.6 sec. 

In addition to TIP and TESOS, we also recorded Ca n K 
3968 A filtergrams with a pixel size of 0."123, a temporal ca- 
dence of 4.93 seconds and a total field of view of 110"x 110". 

The seeing conditions were good to moderate during the 
observations. The ro parameter was around 10, indicating 
spatial resolution no more than 0."5 in the blue range of the 
spectrum at 4500 A, and the resolution around 2" in the in- 
frared. Due to jumps of the Adaptive Optics (AO) system, 
the observed area was lost a few times during the obser- 
vations and we could not completely recover its position. 
Thus, for the further analysis we used only five strictly co- 
spatial TIP scans (from 2nd to 6th), with a total duration 
of 34 min 41 sec. 
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Fig. 1. An overview of the observations. The panels from left to right are: Ban 4554 A continuum intensity in the units 
of its spatially average value; mask applied to locate granules and intergranular lanes; Ban 4554 A line centre intensity 
in the units of its spatially average value; magnetic field strength from the inversion of Fe i IR lines. Contours (same in 
all panels) mark the locations of granular areas with magnetic field above 1.2 kG. 

In summary, after all the reduction and alignment pro- 
cesses, our observational material consists of the following 
datasets: 

— Five TIP maps of Stokes spectra of Fe i lines at 1.56 /zm, 
separated 6 min 50 sec in time and 5."5xl8."5 in size; 

— Time series of Ban 4554 A line spectra over a 
5. "5x18. "5 area, co-spatial with TIP, with a temporal 
cadence of 25.6 sec and a duration of 34 min 41 sec; 

— Time series of Can K filtergrams over 5."5xl8."5 area, 
co-spatial with TIP and TESOS, with a temporal ca- 
dence of 4.93 sec and the same duration as the TIP and 
TESOS scries. 

Figure [1] gives an example of the co-spatial maps of the 
continuum and line core intensity in the Ban line, the mask 
applied to select granular and intergranular areas, and the 
magnetic field strength derived from Fe i infrared (IR) spec- 
tra (see the details below). In the following we do not use 
Ca n K data, but we retain the complete description of our 
observational dataset for a forthcoming paper. The compar- 
ison of the Ba n line core intensity maps and magnetic field 
strength map shows that areas with larger magnetic field 
are, generally, brighter. The Ban line core intensity traces 
magnetic field concentrations in a similar way as Can K 
line does. 

2.2. Calculation of magnetic field from IR Fei lines 

Stokes parameters of the Fei 15648 and 15652 A 
lines were inverted using the SIR inversion code 
( Stokes Inversion based on Resp once functions, see 
iRuiz Cobo fc del Toro Iniestal . 119921) . We fitted aU four 
Stokes parameters of the two Fe i lines. We assumed a one- 
component model including magnetic field. The free param- 



2.1. Data reduction 

As a first step of the reduction process, the data from all 
three cameras (TIP, TESOS and Can K) were corrected 
for the dark current and flat field, following the standard 
procedure. In addition, spcctropolarimctric data from TIP 
were corrected for the instrumental polarization, calibrated 
and de-modulated to give Stokes I, Q, U and V profiles, 
using the standard reduction routines for this instrument. 

We then performed the spatial alignment between all 
three data channels. Here we should recall that TIP records 
all Stokes spectra simultaneously, however each spatial slit 
position is separated 27.3 sec in time. In contrast, TESOS 
records all spatial points simultaneously, but each spectral 
position is separated 25.6/37 sec in time. The alignment of 
these data was far from trivial. 

Firstly, we compared the images of a reticle target, 
scanned during the observations by both instruments. From 
them we calculated the rotation angle and scaling factor 
to apply to TESOS data in order to have the same ori- 
entation and spatial sampling as TIP. We then rotated 
and re-sampled TESOS data accordingly. Secondly, we con- 
structed "false" continuum Stokes / maps from TIP (keep- 
ing in mind that they are not strictly simultaneous in time) , 
and continuum Ban images from TESOS. To make a better 
match, we artificially decreased the spatial resolution of the 
TESOS continuum maps to make it similar to the resolution 
of the TIP maps in the infrared. The TIP continuum image 
was "moved" over the TESOS continuum image to find the 
location with the best correlation. This has allowed us to 
align both images with a 1-pixel precision of TIP (0."185). 
The common area between both instruments had the size of 
5."5xl8."5. A similar procedure was applied to the Can K 
data, using Ban line core images as a reference to find the 
best correlation. 
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Fig. 2. Histograms of the magnetic field strength from the inversion of the Fel hnes. Left: histogram over ah spatial 
points; middle: only pixels with Ban continuum contrast larger than the average (granules); right: only pixels with Ban 
continuum contrast lower than the average (intergranular lanes). 



eters of the inversion were: temperature (with a maximum 
of 5 nodes), line of sight and microturbulent velocities and 
magnetic field. The latter was assumed vertical and con- 
stant with height. The contribution of a non- magnetic sur- 
rounding was taken into account by varying the stray light 
fraction. An average Stokes / profile over the non-magnetic 
areas was taken as a stray light profile. The rightmost panel 
of Fig. [T] gives the magnetic field strength of the first (of 
the five) TIP maps used in the current study. The observed 
facular region contained features with a rather strong mag- 
netic field. In fact, the magnetic field strength measured 
directly from the splitting of the infrared Fe i lines was very 
close to the one obtained from the inversion. The average 
stray light factor reached about 30%, and the average flux 
over the area made about 300 G. The histogram of the mag- 
netic field strength measured in all spatial points in the five 
TIP scans is given in Fig. [2] (left panel). While there are 
many pixels with intermediate field strength (300-800 G), 
the histogram peaks at 1300 G. 



2.3. Calculation of velocity and intensity variations from Ball 

The velocity and intensity variations were calculated from 
the Ban line profiles followin g the A-meter technique 
see for the detailed description Stebbins fc Good^ . Il987t 



Shchukina et al 



Shchukina et al 



[20091: iKostik et al.Ll2009D ." As follows from 
()2009[ ). the formation height of Ban line 



spans all the photosphere from to almost 700 km. 

We used 14 reference widths of the average Ban 4554 A 
line intensity profile over the whole time and space. For this 
average profile, and for prescribed line widths, W , we cal- 
culated the intensity value, I{W)., and_the displacements 
(velocities) in the red and blue wing, yr(W^) and Vb{W), 
corresponding to each width W . We repeated the same pro- 
cedure for each individual Ban profile at each spatial posi- 
tion and time moment t. We then found intensity variations 
5I{t,x,W) and red and blue wing velocities Vr{t,x,W), 
Vb[t,x,W) at the corresponding 14 reference widths. We 
are interested in fluctuations of these quantities given by 
the following relations: 



5I{t,x,W) 

5Vr{t,X,W) 

5Vb{t,x,W) 



I{t,x,W) -I[W) 

Vr{t,X,W)-Vr{W) 

Vb{t,x,W)-Vi,{W) 



(1) 



The intensity variations are normalized to their correspond- 
ing mean value for each W level, to give the contrast vari- 
ations, (5C, as: 



5C{t, X, W) = lit, X, W) /i(W) - 1 



(2) 



We would like to remind here that, applying the A-meter 
technique, both intensity (contrast) and velocity fluctuate 
for a given fixed reference width W , since this width can 
correspond to a higher or lower section of the profile, and 
can be displaced in wavelength due to velocities. 

The zero velocity reference was placed requesting 
{5Vr{t,x,W)) = and {5Vb{t,x,W)) = 0, averaged over 
time and space for the reference width W closest to 
the line core. This determination of zero velocity refer- 
ence is less affected by the standard convective blue shift 
than a granulation-averaged intensity profile, because to 
first order the granular intensity-velocity correlation does 
not infiuence the averaged line minimum, as long as the 
granulation is resolved. According to several studies (e.g. 
iLopresto fc Pierc^ |1985() the lines forming high at the at- 
mosphere (near the temperature minimum, as our Ban 
line) have a very little convective blue shift. The above cited 
work does not provide information about the Ba ii 4554 A 
line, however for lines formed at similar heights (Fe i 3767 A, 
Nai 5896 A, Ki 7699 A), the remaining blueshift is around 
60 m s~^. Give than, we realize that our determination of 
the zero velicity position may contain a blue shft, but its 
value is not expected to be larger than 60 m s~^. 

We should keep in mind that the formation heights of 
the blue and red wing intensity points of the line profile 
are, in gene ral, not the same because of the Doppler velocity 
shift (iShchukina et al.l , l2009l ). so the red and blue velocities 
may correspond to different heights. Despite this, in the rest 
of the paper we use average velocities: 



5V{t,x,W) = {5Vr{t,x,W) + 5Vb{t,x,W)) /2 



(3) 



The formation heights corresponding to each width level 
are expected to vary significantly in space, being different 
in granules and intergranular lanes. To ascribe an approx- 
imate formation height to each of the lA W levels, we fol- 
lowed NLTE (non-Local Thermodynamic Equilibrium) ra- 
diative transfer c alculatio ns, described in previous works 
(jShchukina et all . .2009; Kostik et al.l . I2009D . In the follow- 
ing, we will refer to "continuum" velocity and "continuum" 
intensity (contrast) in Ban as velocities and intensities 
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corresponding to the highest part of the profile, i.e., the 
deepest layer in the atmosphere sampled by our procedure. 
Through the paper we adopt a positive sign for upflowing 
velocity (motion toward the observer). 

2.4. Separation of convective and wave components 

Variations of intensity 61{t, x, W) and velocities 
SV{t, X, W) are mainly caused by oscillatory and convective 
motions. Since here we are primary interested in convec- 
tion, we filtered out oscillatory motions in the Fourier 
space based on the k — uj d i agram (iKhomenko et al.l . 120011 : 
iKostik fc KhomenkoL 120071: iKostik et all . l2009t) . 



2.5. Separation into granules and intergranular lanes 

We fragmented all spatial points into granular and inter- 
granular points using as a criterion their continuum in- 
tensity in Ban line. We removed from the consideration 
the dark pores (located mostly in intergranular lanes). 
Intergranular bright points, characteristic for facular areas 
at high spatial resolution (e.g. at 0."1), are not distinguish- 
able at the resolution of our observations (see Fig. [1]). The 
second right panel of this figure gives an example of the 
applied granular mask. 



3. Results 

To study the relationship between the magnetic field, veloc- 
ity and intensity in our data we proceeded in the following 
way. The magnetic field values, measured in all pixels and 
in all five maps, were sorted by their magnitude, and then 
divided into 20 groups (or bins), each of them containing 
an equal number of pixels. Each bin was assigned a value of 
the magnetic field strength by averaging the B values over 
all points it contains. The Ban intensity and velocity values 
(at different levels of the profile) were measured and aver- 
aged separately for each group. In addition to that, we have 
performed a similar procedure, but separately for granu- 
lar and intergranular points, previously segmented from all 
data. 

Figure [3] (left panel) shows values of the unsigned con- 
vective velocities observed in Ban line at different heights 
in the solar atmosphere, binned for intervals of varying 
magnetic field strength. The dependence on the magnetic 
field strength is not monotonic. The velocity initially de- 
creases with increasing magnetic field strength, reaching a 
minimum near B = 1000 G, and then increases again. A 
similar behaviour of the velocities is observed at all of 14 
levels of the Ban line profile. For a more clear view, the 
figure only presents data for three levels, corresponding, 
approximately, to heights of 0, 400, and 650 km. Note that 
unsigned velocities are represented and, thus, such a be- 
haviour is not surprising. It can be easily understood in 
terms of magneto-convection, taking as example the be- 
haviour of the quiet Su n internetwork magnetic fields (e.g., 
IKhomenko et al.l . l2003l see their Figure 7). Indeed, as our 
observations are taken close to the disk centre, the unsigned 
(absolute) values of the vertical velocities are expected to 
be larger in granules (with weaker field) and in intergran- 
ular lanes (generally containing strong magnetic field con- 
centrations), being small at the transition between granules 
and lanes (where the field is of intermediate strength). The 
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Fig. 4. Convective contrast measured from the Ban line 
at the deepest atmospheric layer as a function of magnetic 
field strength from Fe i lines. Black dots: all pixels are con- 
sidered; stars: only granular points; open circles: only in- 
tergranular points. 



existence of regions with almost zero component of the ver- 
tical velocity at the transitions between granules and lanes 
was recently reported from the quiet Sun data taken b y 
IM aX/SuNRiSE inst r umeri t (iMartmez Fillet et all 1201 If ). 
see IKhomenko et al.l (|2010t ). 

The middle and right panels of Fig. [3] show the be- 
haviour of the Ban velocities separately for granules and 
lanes, this time only for the deepest "continuum" atmo- 
spheric layer and taking into account their sign. As can be 
seen, upward velocities of granules decrease with increas- 
ing magnetic field strength, and the downward velocities in 
intergranular lanes are increasing, in agreement with what 
was said above. 

Figure |4] shows the contrast variation measured at the 
deepest "continuum" atmospheric layer from the Ba ii pro- 
file binned for magnetic field intervals of varying strength. 
The contrast is defined as a standard deviation of the fluc- 
tuations of the quantity SC, see Equation Here we com- 
pare three different statistics: (1) all pixels are considered 
(black dots); (2) only granular points are considered (stars); 
and (3) only intergranular points are considered (open cir- 
cles). The contrast is given in per cents with respect to the 
mean value over the whole observed area and time, so that 
dark areas have negative contrast and bright areas have 
positive contrast. In general, in all the cases, the contrast 
decreases with magnetic field strength. When considering 
all the points, the contrast changes from positive to nega- 
tive as the magnetic field strength increases. The reason is 
similar to the case of velocities, i.e. brighter granules have 
generally weaker field than darker intergranular lanes, same 
as for the magnetic fields in the quiet Sun. At the same 
time, it is interesting that the dependence on contrast be- 
comes much weaker if we consider statistics separately for 
granules and lanes. For lanes, the dependence is within the 
error bar limit, meaning that the contrast of intergranular 
lanes is nearly independent of the field strength. We will 
discuss the implications of this finding later in Section H) 

Another known feature of solar granulation in the quiet 
non-magnetic areas is its contrast sign reversal and the 
velocity direction reversal observed at the middle photo- 
sphere. Here we investigate this feature, but for the facular 
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Fig. 3. Convective velocities measured in Ball line as a function of magnetic field strength obtained from the Fei lines. 
Left panel: unsigned values of the velocity, all pixels are considered. The three group of points correspond to the velocities 
at three levels of the Ban profile, larger points mean higher in the atmosphere (0, 400 and 650 km, approximately). 
Middle: only granular points, velocities at the deepest level of the atmosphere. Right panel: only intcrgranular points. 
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Fig. 5. Upper panels: height of the velocity sign reversal as a function of "continuum" velocity (left) and "continuum" 
contrast (right) from the Bail profile. The mean value is 330±160 km. Each symbol is an average over a bin with an 
equal number of data points. Bottom panels: height of the contrast sign reversal as a function of the "continuum" velocity 
(left) and "continuum" contrast (right) from the Ban profile. The average value is 420±130 km. 



area, searching for dependences of the heights of reversal 
on the strength of gran ular motions at the continuum for- 
mation level, similar to iKostik et al.l (j2009t) . 

In our data, only about 37% of elements show the con- 
trast reversal with height, and 26% show the reversal of the 
direction of motion, as calculated from the Ban A- meter 
velocities and intensities. We only used these elements for 
the results presented in Fig. [51 The scatter of heights of 
the velocity and contrast sign reversals is very high. These 
heights vary in a range from 50 to 650 km, i.e. spanning al- 
most the whole formation height range of the Ba n line pro- 
file. Nevertheless, if these heights are plotted as a function 



of "continuum" velocity or "continuum" contrast, there is 
a clear dependence, si milar to the case of the quiet Sun 
Ban line data used in iKostik et al.l (|2009l ). Fig. [5] shows 
this result. We calculate the error bars as a standard de- 
viation over the values for each bin, and then average over 
all bins. The larger is the continuum velocity and contrast 
of convective elements, the larger is the height where they 
change their direction of motion (upper two panels). On 
average, this height is about 330 ± 160 km. Same happens 
for the contrast sign reversal (bottom two panels) . The de- 
pendence on the continuum contrast is more pronounced 
than on the velocity. On average, the contrast sign reversal 
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Fig. 6. Heights of the velocity and contrast sign reversal as 
a function of magnetic field strength. 



of the convective elements takes place at about 420 ±130 
km. Interestingly, both average reversal heig hts are about 
100 km larger than in the quiet Sun data of iKostik et al.l 
(|2009D («210 and wSSO km, compared to 330 and 420 km). 

Both th e stat istics over the quiet Sun data of 
iKostik et al.l ()2009() . and the facular data in the present 
study suggest that the velocity sign reversal happens at 
lower heights that the contrast sign reversal. This should 
not be interpreted in a way that the same ascending element 
first changes the direction of motion at one height h\ and 
then the contrast at another height /12, where h\ < /12. We 
do the statistics for the contrast and velocity sign reversals 
independently. Thus, we can have independent elements, 
one changing the direction of motion at h\ and another 
one changing the contrast at /12. Also, the descending ele- 
ments can first change their contrast and then the direction 
of motion at a lower height. We do not differentiate between 
the descending and ascending elements, neither between 
the bright and dark elements when calculating the reversal 
heights. Note also, that only about small percentage of ele- 
ments reverse their contrast and direction of motion in our 
data. 

Our next step was to investigate if the reversal height 
depends on the magnetic field strength. To increase the 
statistics, we averaged indistinguishably velocity and con- 
trast reversal heights in all the data points for bins of vary- 
ing magnetic field strength. The result is given in Figure El 
As can be seen from the figure, there is a hint for such a 
dependence within the error bar limit. The reversal height 
is roughly independent of the magnetic field strength till 
about 1000 G. For the larger field strengths, the reversal 
height increases from about 350 to 430 km. 

Yet another important parameter characterizing gran- 
ulation is the correlation between intensity (contrast) and 
velocity of granular elements. This correlation is expected 
to be maximum at the continuum formation layer, de- 
creasing to zero in the middle photosphere (see e.g., 
iKostvk fc Shchukinal . l2004l and references therein). The 
values of the correlation measured by the latter authors do 
not exceed 75% even in the deepest layer, meaning that 
total coincidence of the intensity and velocity maxima is 
rare for the solar photosphere. Here we investigate how the 
correlation coefficient between the velocity and intensity 
changes as a function of the magnetic field. Figure [7] shows 
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Fig. 7. Top: correlation coefficient between the "contin- 
uum" intensity and "continuum" velocity measured from 
the Ban line profiles as a function of the magnetic field 
strength. Bottom: same for the velocities and intensities at 
the highest level ~650 km. 

this dependence for the two atmospheric levels: the deep- 
est "continuum" levels (top) and the highest level around 
the Ban line core (~650 km, bottom). It appears that at 
the "continuum" level, the correlation is maximum for the 
intermediate values of the magnetic field strength (500-700 
G), reaching as high as 77%. The correlation coefficient de- 
creases down to 60% as the magnetic field gets stronger. 
The variations between the maximum and minimum values 
are larger than the statistical error bar (also shown in the 
Figure). Unlike that, the correlation at 650 km is negative 
for intermediate magnetic field strength, characteristic for 
inverse granulation. For larger field strengths above 1000 G, 
the correlation turns positive reaching 20% at maximum. 
This result indicates that strong magnetic field in facular 
areas help preserving the structure of convective elements 
till larger heights. 

4. Discussion 

Before discussing our results, let us recall the medium 
spatial resolution ('^0."5 in the blue part of the spec- 
trum) of our observational data which, in principle, could 
affect our results. In particular, we do not "see" man- 
ifestations of strong kG magnetic field concentrations 
expected to be a ppear as b r ight p oints (BPs) in in- 
tergranu lar lanes (Title et al.l . 119921: iBerger et atl . 120041: 
Narav an &: Scharm er. 2010). Thus, there might be a chance 
that after segmenting all data points into granular and in- 
tergranular points, unresolved BPs may appear belonging 
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Fig. 8. Convective motions in the solar atmosphere accord- 
ing to 16-column modeL Upper panel shows the number of 
cases corresponding to each of the 16 types of motions (as 
defined at the bottom panel) between and 650 km. 



to either as granular or intergranular groups, depending on 
the atmospheric blurring. Possibly, the second maximum (B 
= 1300 G) on the histogram of the magnetic field strength 
distribution over granules (Fig. [2j middle panel) is due to 
unresolved such BPs. For the purpose of visual compari- 
son the granular pixels with magnetic field above 1.2 kG 
are marked red in Fig. [T] At least visually, there is nothing 
particular about the location of such pixels. To check if the 
unresolved BPs affect our results, we have repeated all the 
relevant calculations presented above, but this time exclud- 
ing from consideration the granular points with magnetic 
field strength above 1.1 kG (a value corresponding to a lo- 
cal minimum on the histogram in Fig. [2j middle panel) . We 
have not found any significant differences in our results. 

The dependence of convective velocities on magnetic 
field strength can be compared to the results of simi- 
lar studies by other authors. Measurements for the plage 
and facular areas are not that frequent and the most 
adequate compa, ri son can be done with the works of 
iTitle et al.l (ll992l):lMontagne et al.l (ll996D : lMo"rinaga et al.l 
(j2008[) : iNaravanfc Scharmeij (j2010D . These authors used 
different method for obtaining magnetic field strength from 
their measurements (inversion of spcctropolarimetric data 
or magnetogram signals). Also, data of different spatial res- 
olution (between 0."15 and l."0) and at different wave- 
lengths (spectral lines from 4305 A to 6768 A) are pre- 
sented. Thus direct one-to-one comparison with our results 
is not possible. Nevertheless, we can state that there is a 
qualitative agreement between our results and those from 
the works mentioned above, as for velocity distributions. 
The granular motions in plage areas on average exhibit up- 
flows for intermediate field strength below 600 G. For the 
larger fields above 600 G the structures show tendency for 
downfiows (except for iTitle et all Il992l where downflows 
are dominant for the magnetogram signal above 100 G). 
Our res ults presented in Fig, p better agree with the re- 
sults by iMorinaga et al.l (|2008| )7 though the latter authors 



study the velocity flows in the vicinity of a solar pore, unlike 

the quiet plage analyz ed here. 

Evans fc Catalanol ()1972D . and iHolweger fc Kneeil 
(1 19891 ). were apoarentlv the first ones to discover the 
contrast sign reversal with height in the sola r atm osphere 



Later, thi s phenomenon w as inves tig ated bv [ Kuccra et al 



fl99i; 
(.2004): 



Espagnct e t al.l (I1995D: iKostvk fc Shchukina 
Puschmann et al.l (I2005D: IStodilkal J2006D : 



iJanssen fc Cauzzil ^2oMh ICheung et all (|2007l ) have 



shown from numerical MHD simulations, that the tem- 
perature (and, thus contrast) reversal in the photosphere 
is due to the competition between expansion and cooling 
of rising fluid elements and their radiative heating. The 
heights of contrast reversal derived from observations by 
different authors vary in a wide range from 60 to 350 km. 
Kostik ct al. (200^ pointed out that this scatter may be 
attributed to the difference in the formation heights of 
the spectral lines used in different observational studies. 
More recently, it was found that not only contrast reverses 
sign in the solar atmosphere, but also the direction 
of motion of convec ti ve el e ments revers e s with height 
(IKostvk fc ShchukinaL 120041: iKostik et al.1 . I2009D . These 
studies were performed for granulation observed in very 
quiet non-magnetic regions. Our results from Fig. [5] point 
out that both contrast and velocity sign reversal take place 
in the facular regions as well. The moderate magnetic fleld 
in facula (400-900 G), does not affect the height where the 
reversal takes place compared to the quiet solar regions 
(Fig. [6]). However, stronger field above 1 kG "forces" the 
convective elements to travel over larger vertical distances 
without changing the sign of contrast and direction of 
motion. Apparently, strong magnetic field in facular areas 
promotes more effective energy transfer in the upper layers 
of the solar atmosphere by convective elements, since they 
can reach larger heights. At the first instance, it might 
seem that this conclusion is in disagreement with the 
results shown in the top panel of Figure |7| regarding the 
correlation between velocity and contrast. The correla- 
tion coefficient almost linearly decreases with increasing 
magnetic field. However, the correlation at the heights of 
formation of the Ban line core behaves in the opposite 
way (bottom panel of Figure [7]). While the absolute values 
of the correlation in the upper atmosphere are lower, as 
expected, the correlation is, on average, positive for larger 
magnetic field strengths above 1 kG, indicating that the 
structure of convective elements is maintained till that 
height. 

Yet another prove of the latter result is provided by 
Figure m This figure compar es the types of moti on observed 
in the quiet area (data from lKostik et al.l . l2009D and m our 
facular area. The bottom panel of the figure indicates our 
classification of the types of motion: sign indicates con- 
vective elements moving upwards or those whose contrast 
is above the average; sign "-" indicated downward moving 
elements or those with the contrast below the average. We 
compare the velocity sign and contrast of convective ele- 
ments at the "continuum" level and the at highest level 
available for us around the formation of the Ban line core. 
Interestingly, we find that the number of convective ele- 
ments that do not change neither contrast nor velocity di- 
rection in the facular area to be about 53% (motions of type 
1 to 4). Compared to this, in the quiet area only 31% of 
convective elements behave this way. We conclude that the 
strong magnetic field of the facular region helps to stabilize 
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the convective elements and allow the convective energy to 
reach larger heights. 

Possibly the most interesting result of our study is 
the one presented in Fig. H) The general decrease of the 
contrast with increasing magnetic field (if all the points 
are considered indistinguishably) seems to be a natu 
ral result and confir r ns the resul t s by the earlier stud- 
ies (iTitle et al.l. 119921: iTopka et all. Il992t iMontaene et al 



1996 1 iTopka et al.l . 119971: iBerger et al.l l2007l : iKobel et al 
20111) . However, the behaviour of intergranular pixels, con- 
sidered separately from the rest, is significantly differ- 
ent (dashed line in Figure 2]). Namely, the intergranular 
contrast is essentially independent of the magnetic field 
strength within the err or bar limit. The sa me conclusion 
was reached earlier by iBerger et al.l (|2007[ ). The authors 
claimed that in facula regions close to the limb one rather 
detects granular walls, that should be essentially indepen- 
dent of the magnetic field strength. Here we confirm this 
finding but for the facular region close to the disk centre. 

The decrease of the contrast with increasing magnetic 
field obtained in previous studies seems to be a result of 
the particular distribution of the magnetic field strength 
in these structures (Fig. [2]). The granules have generally 
weaker fields, while intergranules have stronger ones. By 
averaging indistinguishably the intensity values for differ- 
ent B intervals, one would naturally obtain a larger num- 
ber of granular pixels for the bins with lower B (i.e., larger 
contrast); and the opposite for intergranules. Thus, it is im- 
portant to segment granular and intergranular pixels prior 
to analyzing the dependence of the contrast on B. 

All the statistical dependences found in our paper are 
more significant that the corresponding error bars, given in 
all the figures. The only exception is the dependence shown 
in Fig.[Bl which is within the error bar. The error bars were 
calculated as standard deviation over the values averaged 
for a given bin. 

The advantage of our work over the similar studies is 
the real magnetic field measurements obtained from the in- 
version of extremely magnetically sensitive Fe i spectra in 
the infrared, instead of using a magnetogram signal. The 
alignment of our data sets from both instruments was done 
with a 1 pixel precision (0."185), improving the statistical 
significance of our results. It is important to recall, how- 
ever, that the diffraction limits in the visible and in the 
infrared are very different, leading to a naturally lower spa- 
tial resolution of the infrared data. For a better alignment 
we artificially reduced the spatial resolution of the contin- 
uum Ba II maps to make it comparable with the one of TIP 
Stokes / maps. The lower resolution of TIP data might, 
in principle affect the precise value of the magnetic field 
derived from this data, since the magnetic field patches 
are more diffused and magnetic elements of mixed polari- 
ties might cancel. However, the cancelations are expected 
to be more important in the mixed-polarity inter-network 
areas, and not in a (mostly unipolar) plage area as the 
observed one. Since we derive the magnetic field from the 
inversion of Stokes profiles, using a stray-light component 
as a free parameter, this improves our accuracy in the mag- 
netic field strength determination. The magnetic field val- 
ues given at the horizontal axis of Figs. El SI IHl and [7] are 
averages over the intensity bins. While these average val- 
ues might, in principle, be affected by the cancelations in 
a mixed-polarity area, we do not expect this to be impor- 
tant in a uni-polar plage. Another possible source of error 



is that our data do not have enough resolution to resolve 
intergranular bright points. However, as pointed out at the 
beginning of this section, we have cross-checked that re- 
moving the "granular" points with _B > 1.2 kG from our 
analysis (as possible unresolved bright points) provides the 
same results, in particular the one given in Fig. S) This 
confirms once more our conclusion that the contrast of in- 
tergranular lan es is independent on the magnetic field, in 
agreement with IBerger et al.l ()2007l ). 

It will be interesting in the future to check the results 
of our work using data with a larger spatial resolution and 
maintaining the advantages of the simultaneous Stokes po- 
larimetry and velocity /intensity measurements till the up- 
per photosphere. 



5. Conclusions 

Using simultaneous spectropolarimetric observations in the 
infrared Fei lines at 1.56 fim and spectral observations of 
the strong line of Ban 4554 A formed high in the atmo- 
sphere, with a moderate temporal resolution, we have inves- 
tigated the behaviour of the convective component of the 
intensity and velocity fields in the facular region close to 
the solar disk centr. The Ban line allowed us to perform a 
study over a wide range of heights, from the continuum for- 
mation height up to about 650 km in the solar atmosphere. 
We have reached the following conclusions: 

1. At the bottom photosphere, the convective velocities of 
granules decrease with magnetic field strength, while 
the convective velocities of intergranules increase with 
magnetic field strength. 

2. Similar to the quiet regions, we detect that in facular 
regions the contrast of granulation and the velocities 
of convective elements reverse their sign with height. 
The height where such reversal takes place depends 
on the strength of convective elements at the bottom 
photosphere, but also on the magnetic field strength. 
Convective elements above stronger magnetic features 
reach higher in the atmosphere without breaking. 

3. The correlation coefficient between velocity and inten- 
sity at the bottom-most level decreases with the mag- 
netic field. In the upper atmosphere, above 500 km, 
the correlation gets larger with increasing the magnetic 
field. 

4. The contrast of convective elements decreases with in- 
creasing magnetic field strength. But if intergranular 
lanes are considered separately, their contrast is nearly 
independent of the field strength. 
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